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Abstract . Shock wave data, thermodynamic and phase diagram data for ice, 
porous ice, and water are taken together with a Rice-Walsh-Bakanova equation of. 
state to define the shock pressures and impact velocities required to induce 
incipient melting (IM) (6 GPa), complete melting (CM) (10 GPa), and passage 
through the vapor-liquid critical point (CP) upon isentropic release (22.5 GPa). 
Upon expanding along the isentrope which passes through CP ~0.5 mass fraction 
vaporization occurs until a pressure of ~61 kPa (61 mbar) is achieved. Below 
this pressure ice sublimates and ~0.4 mass fraction H 2 0 gas is in equilibrium 
with ice I. The minimum impact velocity required to induce IM, CM, and isentro- 
pic release through CP is 2.1, 3.0, and 4.5 km/sec for silicate impactors. For icy 
projectiles Hugoniot states depend only weakly on initial temperature of ice. The 
IM, CM, and CP isentropes are achieved upon impacting an icy surface at veloci- 
ties 3.4, 4.4, and 7.2 km/sec, respectively. We observe that a partial H 2 0 pres- 
sure below 61 kPa and temperatures below 273K, ice partially vaporizes and 
requires 2800 kJ/kg, of heat of vaporization for complete sublimation. We exam- 
ine the hypothesis that the smaller satellites of Saturn having mean densities in 
the 1.1 to 1.4 Mg/m 3 range represent primordial accreted planetesimal conden- 
sates (60% (wt.) H 2 0, and 40% (wt.) silicate) formed in the proto-Jovian and 
Saturnian accretionary planetary discs. These densities are in the range 
expected for water-ice/silicate mixtures constrained to the solar values of 0/Si 
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and 0/Mg atomic ratios. We demonstrate that if the large satellites accreted 
from the same group of planetesimals which formed the small Saturnian satel- 
lites impact vaporization of water upon accretion in a porous regolith, at low H 2 0 
partial pressure, can account for the increase in mean planetesimal density 
from 1.6 Mg/m 3 (43% H 2 0 + 57% silicate) to a mean planetary density of 1.9 
Mg/m 3 for Ganymedean-sized water-silicate objects. If impact volatilization of 
initially porous planetesimals is assumed, we demonstrate that starting with 
planetesimals composed of 54% H 2 0 and 40% silicate (1.35 Mg/m 3 ) partial devola- 
tilization upon accretion will yield a Ganymede-sized planet, having a radius of 
2600 km and a density of 1.85 kg/m 3 , similar to that of Ganymede, Callisto, and 
Titan. 

1. INTRODUCTION 

We apply knowledge of the properties of water and ice over a wide range of pres- 
sures and temperatures and describe constraints on the shock vaporization 
processes for water and ice in the solar system. In particular, we examine the 
role of impact vaporization acting during the formation of the Jovian and Satur- 
nian satellites in an attempt to explain the observed density in terms of compo- 
sition of these rock and ice objects. The imaging instruments on the Voyager I 
and II spacecraft (Smith et al., 1979a, b, 1981, Morrison, 1982) have provided 
definitive knowledge of the volumes, and hence, densities, as well as, the surface 
characteristics of at least 12 satellites of Jupiter and Saturn having planetary 
radii ~ 100 km (Table 1). 

Pollack and Reynolds (1974), and more recently, Lunine and Stevenson 
(1982), and Pollack and Finale (1982) have discussed the origin of the larger 
Jovian satellites within a framework of condensation from a "protostellar" disc of 
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gaseous and particulate material rotating around a newly formed proto-Jupiter 
(of approximately solar composition). This disc is much like the larger scale 
protosolar disc from which the sun and the planets themselves form. Like the 
solar system, itself, this type of model of the Jovian satellite system implies a 
higher temperature close to Jupiter and thus ultimately less relative condensa- 
tion of volatiles (e.g. HaO, NH 3 ) in the satellites close to the planet as compared 
to the outer satellites. The Jovian system clearly suggests such a formational 
history, as the inner satellites (Amalthea, Io, and Europa) have higher mean den- 
sities, 3.0 to 3.6 Mg/m 3 , and in an analogy to the terrestrial planets) contain a 
higher fraction of silicates than do the outer larger (lower density) satellites 
(Ganymede and Callisto). The latter have densities of 1.8 to 1.9 Mg/m 3 (Fig. 1). 
In such formation models of the Jovian satellites, it is usually assumed that 
these accreted from small solid objects which formed via gravitational or other 
instabilities in the proto-Jovian disc. The inner, satellites accumulated lower 
inventories of volatiles as a result of either or both of the following two effects: 

(1) Because of the higher temperatures in the vicinity of the early Jupiter, 
relatively less volatiles can condense. For example, in the Lunine and Stevenson 
model the temperatures and pressures at the orbits of Amalthea and Callisto 
are 1350K and 13 MPa and 150K and 0.01 MPa, and 

(2) Because the Keplerian velocity around Jupiter of the proto-satellites 
and hence the speed with which they overtake and accrete smaller solid objects 
imbedded in the gaseous nebula have a lower encounter velocity with increasing 
radius from Jupiter, less impact vaporization can occur during accretion. 

Whereas the "miniature solar system" model for the Jovian satellite system 
which has relatively more volatile material condensing with increasing radius in 
one form or other should be successful in explaining the high density inner 
satellites versus the low density outer satellites. In the case of the Saturnian 
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system, the satellite mean density versus radius from the planet relation is 
irregular. As can be seen in Fig. 1, there is no obvious correlation of mean satel- 
lite density with orbital radius. 

Examination of Table 1 suggests that one could classify the Jovian and 
Saturnian satellites into three groups according to mean density: silicate (sil), 
primordial (p), and partially devolatilized (pd) (e.g. Cole, 1984). By silicate den- 
sity we mean the density of terrestrial planetary mantles 3-3.5 Mg/m 3 . An esti- 
mate of a primordial density can be obtained by assuming the solar 0/Si atomic 
ratio of 15.5 (Ross and Aller, 1976). We assume that this ratio determines the 
abundance of H 2 0 which is in turn controlled by the abundance of oxygen rela- 
tive to silicon and hence the abundance of Si. For example, if we consider SiOg 
as being controlled by the abundance of Si then we infer that an atomic ratio is 
required of 5.16 moles of H a 0 to 1 mole of Si0 2 . This ratio thus implies a density 
of 1.23 Mg/m 3 for a water-ice and quartz mixture. Assuming that the most abun- 
dant anhydrous silicate in the solar system is probably enstatite MgSi0 3 and util- 
izing the solar atomic ratio of 0/Mg of 17.38 and assuming the H a O/MgO ratio 
controls the density of primordial low temperature volatile-silicate condensate 
mixture, a density of 1.08 Mg/m 3 is inferred by the same reasoning. The addi- 
tional small amounts of oxidized iron will, of course, increase this density. We 
assume that densities in the above range, correspond to a silicate plus some oxi- 
dized iron plus volatile primordial condensate and this is assumed to be 
representative of densities of the group of the smaller satellites in the Saturnian 
systems indicated by p. 

In the present paper we first examine the shock wave and phase diagram 
data for water, solid and porous water-ice and define the impact velocity 
required for shock-induced melting and vaporization of water and ice at various 


temperatures. 
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We then examine a possible model of accretion of icy satellites which 
predicts that the amount of ice devolatilization is related to planetary size. The 
Lunine and Stevenson (1982) model, predicts final satellite density on the basis 
of differing environment for accretion as a function of satellite radius from 
Jupiter. We assume all the satellites in the Jovian and Saturnian system (except 
J5, Jl, and J2) accreted from the same composition of planetesimal material and 
that the least affected samples of these planetesimals are the smaller satellites 
of Saturn and possibly also the ring materials of both planets. The partially 
devolatilized satellites whose density we attempt to predict (Ganymede, Callisto, 
and Titan) are then assumed to be partially devolatilized, indicated as p.d. in 
Table 1. 

2. Shock Wave Equation of State and Thermodynamic Data for Water, Water-Ice, 
and Porous Ice. 

We need to combine knowledge of the Hugoniot curves for ice targets and 
the calculated release isentrope curves to determine the impact velocities 
required to generate sufficient entropy to melt and vaporize ice. The train of 
logic is as follows: 

(a) Given the measured Hugoniot shock velocity, U 3 , versus particle velo- 
city, Up, data, shock pressure versus shock particle velocity curves may be con- 
structed. These allow prediction via the impedance match method (McQueen et 
al., 1970) of the peak shock pressures induced upon impact of a specific projec- 
tiles (e.g. rock or ice) against solid or porous ice or water targets. The amount 
of material subjected to this peak pressure is in the order of 10° to 10 1 times the 
projectile mass (O'Keefe and Ahrens, 1977). 

(b) The single shock data must then be combined with data obtained by 
double shocking water and/or porous ice to the same pressure but at lesser and 



greater temperatures (and hence lesser and greater specific volumes) as well as 
thermodynamic data to calculate the isentropes in pressure-volume tempera- 
ture space. The tactic we employed is that of Ahrens and O’Keefe (1972). We 
utilize pressure volume-entropy-enthalpy data at a single moderate pressure (l 
GPa) and increasing (high) temperatures to define a series of states with 
increasing entropy to specify the foot of isentropes. 

(c) These isentropes are then projected to higher pressures using a com- 
plete equation of state formalism like that proposed by Bakanova et al. (1976). 
Once the intersection with the water and ice Hugoniots are found we can infer 
from the intersection of the various isentropes with the Hugoniot the entropy 
along the Hugoniot and hence determine from a given shock state whether the 
entropy density is sufficiently high such that upon pressure release shock 
induced melting of vaporization will occur. Of special interest are the isentropes 
that cross the Hugoniot which define the shock pressures required to just pro- 
duce incipient melting (IM), complete melting (CM), and bring water to the 
liquid-vapor (steam) critical point (CP). 

(d) When the calculated release isentropes are taken together with the low 
pressure thermodynamic data for the water-ice steam system, as for example, 
summarized in the form of a Mollier diagram (Bosnjakovic et al., 1970) the mass 
fraction of vapor produced by a given entropy density as a function of partial 
pressure of water vapor in the planetary environment may be inferred. 

(1) Shock wave data. 

Because ice covers the surface of more than 20 of the Jovian and Satur- 
nian satellites (except for J 5 , J t , and S 6 ) and is a major constituent of the inte- 
rior of these objects (and it of course covers 70% of the earth) it is important to 
determine the projectile impact velocities and hence shock pressures required 
to induce melting and/or vaporization. 



-7- 


The available Kugoniot data for water-ice and porous-ice in terms of shock 
velocity, U, particle velocity, u, relations of the form 

U = C 0 + X( u + X 2 uS (l) 

is given in Table 2. Shock pressure and shock density (or specific volume, V) can 
be obtained from the constants in Table 2 via the Rankine Hugoniot shock equa- 
tions for pressure, P, density, p, and internal energy, E, or enthalpy, K, respec- 
tively. 


P-Pq = Po u U (2) 

p = p Q U/(U-u) or V = V 0 (U-u)/U (3) 

E-E 0 = (Po+P)(V Q0 -V)/ 2 or H-H 0 = (P o +P)(V 00 +V)/ 2 (4) 

In equations 2-4, P 0 , p 0 , V Q , E„, and H 0 are the crystal density, specific cry- 
stal volume, internal energy and enthalpy in the initial state. In the present 
paper, energies and enthalpies (and entropy) are measured relative to condi- 
tions at 1 bar and OK. V Q0 is the pre-shock initial specific volume which in the 
case of porous ice will be greater than V Q . For no porosity V 0 = V cc . In the case 
of liquid water, the shock wave data are described via a set of bilinear equations 
by Mitchell and Nellis (1992) who assume a high pressure regime extending from 
4.4 to S3 GPa and Bakanova et al. (1976) who defined the high pressure regime 
from 9 to 55 GPa. Double (reflected) shock experiments on water have also been 
carried out by Walsh and Rice (1957), Bakanova et al. (1976), and Mitchell and 
Nellis (1992), the latter to peak pressures of 230 GPa. In addition, a large 
number of shock wave data for solid ice initially at temperatures of 259 to 263 K 
are available over a pressure range from 0.6 to 50.3 GPa (Fig. 2). These data are 
also closely fit by a linear shock velocity-particle velocity relation. Below 3 GPa, 
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shock waves in solid, ice demonstrate interesting structure involving intermedi- 
ate shock states due to phase transitions. These are described by Gaffney 
(1984) in this volume. 

When the final shock state data summarized in Table 2 is cast into the 
pressure-particle velocity plane (Figure 3) via equation 2, it can be seen that 
liquid water (on account of its greater density) initially has a slightly higher 
shock impedance (larger value of shock pressure at a fixed particle velocity) 
than ice and that this small impedance contrast reverses (on account of more 
shock heating in ice) above ~95 GPa at a particle velocity of 7.7 km/sec. Porous 
ices with lower initial densities have correspondingly lower shock impedance. 

The representation of Hugoniot curves in the pressure-particle velocity 
plane via the impedance match method (McQueen et al., 1970) allows straight- 
forward determination of projectile-target peak shock pressures from the 
impact velocities. To demonstrate how this is carried out we indicate in Fig. 2 
the escape velocity of water-bearing objects in the solar system along the 

ft 

abscissa (from 0.16 km/sec for Mimas to 11.19 km/sec for the earth). In gen- 
eral. the infall velocity is the minimum velocity for which a projectile can impact 
the planetary surface in the case where no atmosphere exists or the projectile is 
large enough not to be affected by its passage through the atmosphere. This 
minimum velocity then can be increased substantially in the case of the satel- 
lites of the major planets by the following three factors which are listed in order 
of importance. 

(1) The effect of gravity focusing by the mass of the planet. 

(2) The effect of heliocentric velocity of the planet, and 

(3) The effect of the orbital velocity about the planet. 

The effect of gravitational focusing by the high infall velocity (escape velo- 
city of Jupiter and Saturn, 60 and 36 km/sec, respectively) of objects on their 
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planetary satellites is demonstrated with the mean impact velocity and flux con- 
centrations given in Table 3. Table 3 assumes that for the Jovian system objects 
have a mean encounter velocity of 8 km/sec. At Jupiter’s heliocentric distance 
many of these objects are cometary (see Shoemaker and Wolf, 1982). 

The pressure-particle velocity Hugoniot of the impactors is plotted back- 
ward centered at the impact velocity. Thus shown in Fig. 2, a projectile impact- 
ing at 11.19 km/sec having the pressure-particle velocity characteristics of ser- 
pentinite, will induce a shock pressure of 41, 63, 90, and 90.5 GPa in 0.35 and 0.6 
Mg/m 3 porous ice, solid ice, and water, respectively. In contrast to the major 
and terrestrial planets the minimum infall velocities of the Jovian and Saturnian 
satellites are much lower as indicated in Fig. 3. For solid ice infalling at the 
escape velocity of Ganymede (2.74 km/sec) a peak pressure of only 4.2 GPa is 
induced in the projectile and target. 

(2) Equation of state for water and ice. 

It was recognized as early as 1957 when Rice and Walsh (1957) carried out 
the first extensive analysis of single and reflected shock data for water that the 
formulation of a complete thermal equation of state according to the Mie- 
Gruneisen equation — so effective for metals (e.g. Rice et al., 1957) was infeasi- 
ble for application to water and other molecular fluids. Instead Rice and Walsh 
proposed on the basis of a series of reflected shock data in which different 
enthalpy states were achieved at a given pressure that the thermodynamic 
quantity (9H/3V) P depends strongly on P (and is effectively independent of V or 
T at constant P). They then defined a thermodynamic quantity £ (P) which we 
will all the Rice-Walsh parameter, as 



av 

3T 


£(P) = c p /(av/9T)p = 


p 


(5) 
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Eq. (5) is analogous to the Mie-Gruneisen equation of state, where a state at 
a given E can be calculated for a given pressure, P, by comparing it to the inter- 
nal energy (E r ), and pressure (P r ) at the same volume V. Here E r and P r refer to 
the reference thermodynamic path such as the principal Hugoniot. The Mie- 
Gruneisen equation-of-state is written as: 

E = E r + V/ 7 (P-P r ) ( 6 ) 

where 

7 s V (aP/ 3 E) V ( 7 ) 

In analogy with the Mie-Gruneisen parameter, the Rice-Walsh parameter is used 
to relate enthalpy H and volume, V, at the same pressure, P, to that along the 
reference path. The reference path chosen here is the principal Hugoniot of 
water denoted by the subscript 1 . Thus the Rice-Walsh equation is 

H(P, V) = Hj (P) + $(P) [V— V,(P)] (8) 

This treatment satisfactorily describes the available Hugoniot and reflected 
shock data to 25 GPa and is also concordant with static compression data for 
water at high temperature as well as high pressure and relatively low tempera- 
ture on the properties of steam. Upon obtaining additional reflected shock data 
to 120 GPa and porous ice Hugoniot data, Bakanova et al. (1976) utilized an 
equivalent and more convenient relation than equation 9 specifying internal 
energy at constant pressure. The Bakanova et al. parameter is defined as 

77 (P) = P (9V/ 9E) P (9) 

and is related to the Rice-Walsh parameter by applying the relation 


(3H/9V)p = (3E/3V)p + P 


( 10 ) 
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to obtain 

v( p) = p/(£-p) (ID 

From equation 9 it foLlows that 

E(P. V) - Ej(P) = -^-(V-V^P)) (12) 

where Ei and V t are again the reference internal energy and volume at pressure, 
P. 

To obtain an equation for an isentrope we substitute for E t (P) 

E t = E 0 + (V 0 -V) P/ 2 (13) 

using the Rankine Hugoniot energy equation (relative to the internal energy, E 0 
at STP) into equation 12 and solve for V. This yields 

v 

V = V, (1+77/2) -V f PdV-E i0 +E 0 + PV 0 /2 / P (14) 

Vl ° 

where we have also substituted for E 

v 

E = - / PdV + E i0 (15) 

V ic 

Here V io and E io , and T io are the specific volume internal energy and temperature 
state which defines the foot of the isentrope at pressure P j0 . 

An equation for the temperature along the isentrope can be obtained from 
the differential relation 

0 = TdS = T (aS/ 3T) P dT + T ( 3S / dP) T dP (16) 

substituting in the first and second term, respectively, and using the 
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thermodynamic relations 

Cp/T = (8S/9T) P (17a) 

(as/ ap) T = -(av/ aT) P (i7b) 

yields the differential equation 

dT/T = dP(aV/aT) P /Cp (18) 

which upon substituting for £ (P) from equation 8 and integrating yields 

T * =T “ esp ^ ffr =T, ° exp / s (19> 

In order to calculate the P-V-T-E isentropic path from equation 14 and 19 
the function £(P) or 7](P) must be determined. Rice and Walsh examined four 
reflected shock data to 23.5 GPa and fit these to 

log io £ = 1.17943 + 0.030338 P(GPa) (20) 

whereas Bakanova et al. fit both porous ice data and water double shock data 
and fit these to the relation 

77 (P) GPa = 0.07 [ 1— exp(-0.08436P)] + 0.044095P exp (-0.048202 P) (21) 

As can be seen in Fig. 4 the 7j(P) function from both studies are similar. The 
peak in 77 in the 15 to 20 GPa range is attributed to the increase pressure- 
induced ionization which occurs in water and is reflected in the rapid increase in 
electrical conductivity (Mitchell and Nellis, 1982) of shocked liquid water. To 
further demonstrate the adequacy of the Bakanova equation of state in describ- 
ing the thermodynamic properties of water to high pressures, we apply it to cal- 
culating shock temperatures for water. Comparison of calculated and measured 
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shock temperatures is carried out in Appendix A. 

In order to calculate the specific release isentrope P-V paths for H 2 0 other 
than the paths given by Bakanova (1976), Fig. 5, we can utilize both the phase 
diagram of water (Fig. 6) and the Mollier diagram (Fig. 7). The phase diagram of 
water sketched in Fig. 6 demonstrates that upon shocking either liquid water or 
ice above ~20 GPa it achieves states above the vapor-liquid critical point (647K, 
22.1 MPa, 3.1 m 3 /Mg). It should be noted that a well-defined critical shock pres- 
sure marking the onset of incipient vaporization occurs upon shock release can- 
not be given. Shock-induced vaporization depends strictly on the partial pres- 
sure of H 2 0 in the environment into which shocked water or ice expands. If the 
partial pressure of H 2 0 is sufficiently low (below the ice-I vapor or liquid phase 
line) shocked water ice will always be vaporized upon sufficient entropy input. 

The existence of the ice I-liquid- vapor triple point at 273K and 6.108 kPa 
implies that shocked ice need not melt first but can bypass the liquid field and 
sublimate provided the pressure is below the ice I-sublimation line which is given 
explicitly in Dorsey (1940) as 

log 10 P (torr) = - 4 8.2312 log 10 T - 1677.006 (10~ 5 ) T + (22) 

120514 (10 _1 °) T 2 -6.757169 

The Mollier diagram for water, ice, and steam shown in Fig. 7 defines the 
critical isentrope (vertical lines in Fig. 7) which can be defined for incipient 
melting (IM), complete melting (CM), and the isentrope which passes through 
the vapor-liquid critical point (CP). These critical isentropes in the pressure 
volume plane are plotted in Fig. 8 for ice at 70 and 263K, the range of tempera- 
tures appropriate for the surfaces of the Jovian and Saturnian icy satellites. 
Whereas the isentropes corresponding to IM and CM produced little vapor until 
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relatively Low pressures are achieved, the release isentrope passing through the 
CP gives rise to ~0.5 mass fraction vapor upon reiease through the critical point 
and upon release to pressures below 61 kPa some 0.4 mass fraction is vaporized 

(Fig- ?)• 

The critical isentropes of Fig. 8 and Table 5 are also shown in the pressure- 
particle velocity plane for solid ice in Fig. 2. Examination of where the IM, CM, 
and CP points fall on the ice Hugoniot (Fig. 2) demonstrates that the minimum 
infall velocities for ice onto an icy satellite so as to produce melting or vaporiza- 
tion can be specified. Silicate projectiles impacting at greater than 2.1 krn/sec 
or icy projectiles impacting 3.4 km/sec are required to induce IM of ice. CM of 
ice requires silicate projectiles impacting at least 3 km/sec or ice projectiles 
impacting at 4.4 km/sec. To obtain 0.4 to 0.5 mass fraction vapor upon passing 
through the CP requires ~4.5 km/sec for silicate impactors or 7.2 km/sec for 
pure ice impactors. The latter velocities are in approximate agreement with 
those of Smoluchowski (1983). 

Thus we conclude that although copious melting and vaporization of water 
ices on icy satellites will occur at meteoroid and cometary impact velocities on 
the Jovian and Saturnian satellites (Table 3), vaporization processes during 
accretion will not be significant unless the partial pressure of water is below 61 
kPa and the temperature is below 273K during accretion of these objects. We 
explore the consequences of accretion of the icy satellites in this environment in 
the next section. 

3. ACCRETION OF ICY SATELLITES 

We have demonstrated the importance of the existence of low partial pres- 
sure of H g 0 (below 61 kPa or 6.1 millibar) for vaporization to be important in 
impact processes during accretion. For accretion under low water partial 
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pressure conditions and low temperatures, impact vaporization is, however, 
readily accomplished because the liquid stability field of water is entirely 
bypassed. Examination of the Mollier diagram of water (Fig. 7) in this low tem- 
perature region indicates that an enthalpy gain of 2800 kJ/kg are required to 
completely vaporize Ice I at pressures below the vapor-liquid ice 1 triple point. 
Rather than assume that the Jovian satellites accreted in a gaseous nebula in 
which the partial pressure of water was at the saturation point at temperatures 
varying from (150 to 230K) as assumed by Lunine and Stevenson (1982), in the 
present paper we will examine the consequences of creating icy satellites via 
infall of planetesimals in a dry, undersaturated environment. 

We assume as an extreme case once material is vaporized upon impact, it is 
lost from the planet. The more complicated problem of precipitation versus 
thermal loss of H 2 0 and the general thermal state of the accreting planetary 
atmosphere is not addressed. 

We shall first examine the case of accretion of a pure icy planet at low pres- 
sures. If we assume that a porous icy regolith forms in this situation we can 
model the fraction of ice impact vaporized, f, as being equal to: 

f = V, 2 / (2 ■ AHsub,) (23) 

where Vj is planetesimal impact velocity. 

We conclude that it is most likely that at low temperatures and at undersa- 
turated conditions a porous ice regolith some ~10°to 10 2 km thick would be 
retained on accreting icy satellites because: (l) Porous regoliths are known to 
adhere to objects as small as ~10 2 km diameter asteroids in the case of Phobos 
and Deimos. (2) At low temperatures ice retains ~50% porosity until overbur- 
den lithostatic stresses of ~5 MPa are exerted (Anderson and Benson, 1963). 
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At any point in the accretion of a uniform density planetary satellite, p 0 , the 
minimum infall velocity, V inijn is equal to the escape velocity, V esc which is 

V esc = VBtt Gp 0 /3 R (24) 

where G is the gravitational constant. We now assume that 

Vi = V i-min = V esc (25) 

and calculate for a pure water planet how much incident material accreted is 
retained versus planetary mass (Fig. 9). We observe that at planetary masses of 
~10 24 g the infall velocity is so slow (0.5 km/sec) that little impact vaporization 
of H 2 0 occurs. However, as the pure H 2 0 satellite accretes in this environment 
an increasing fraction of the mass delivered is impact vaporized and becomes 
lost to the object. As can be seen in Fig. 9 when the planet mass is ~10 29 g, the 
minimum infall velocity is ~2 km/sec and only ''■'0.44 mass fraction of the 
incident planetesimal material has been retained. For a pure H 2 0 planet as 
massive as Titan or Callisto only ~15% of the infalling mass is not impact vapor- 
ized. 

If now the same calculation is carried out for 0.4 mass fraction of silicate, 
0.6 mass fraction water, planetesimals typical of material condensing from the 
primordial Jovian and Saturnian nebulae, and, we assume a silicate density of 
3.0 Mg/m 3 for the silicate fraction, the water vaporization is less extreme as also 
indicated in Fig. 9. Planetesimals with this composition when accreting to form a 
planet retain ~65% of the incident mass for an object when the object becomes 
massive as Titan or Ganymede. 

We sought to examine the question of whether accretion of a water-silicate 
mixture, initial density 1.27 Mg/m 3 which approximates the mean density of the 
smaller satellites (Mimas, Enseladus, Tethys, Dione, Iapetus and Rhea) results in 
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sufficient impact vaporization to obtain a mean density of >1.8 Mg/m 3 
corresponding to larger satellites Callisto, Titan, and Ganymede. As indicated in 
Fig. 10 accretion of such a mixture results in a mean planetary density of only 
~1.5 Mg/m 3 . Moreover, also shown in Fig. 10, the situation is not especially 
improved if one assumes a 10 3 km silicate core for the larger satellites. In this 
case the increase in the average final density is only slight. The interesting 
minimum density at 1.55 Mg/m 3 is caused by the process of initially reducing 
the mean density of the growing planet with a silicate core by accreting water- 
bearing material. However, as the planet grows, an increasing fraction of the 
water-bearing material contributed is shock-vaporized and lost, and finally at a 
radius of ~2300 km the mean density begins to increase again. We conclude 
from these calculations that the simple process of impact vaporization in a 
regolith and a loss of all the shock volatilized water cannot make planets with 
densities as high as Callisto, Titan, and Ganymede from a protosolar mixture of 
60% (mass) water and 40% silicate which both correspond to a primordial water 
silicate mixture as well as approximates the average density of the smaller 
satellites of Saturn. As indicated in Fig. 10 if one relaxes the constraint on the 
planetary composition slightly, and assume 43% H 2 0 and 57% silicate, it is possi- 
ble to produce objects with exactly the density of Callisto, Titan, and Ganymede 
(which correspond to '"'40% water and 60% silicate) by the impact volatilization 
process outlined above. Finally, if instead of assuming a porous regolith, we 
make a slightly more drastic assumption and. assume porous planetesimals (and 
restrict impact vaporization to the mass of planetesimal material only ) it is pos- 
sible to slightly enhance the impact volatilization process to achieve objects 
close to the radius of Callisto, Titan, and Ganymede with material very close to 
the mean density of the smaller satellites of Saturn. In this case our calcula- 
tions demonstrate that starting with 54% H 2 0 and 46% silicate composition 



(starting density 1.35 Mg/m 3 ) planets having the Ganymede radius of 2600 km 
and a density of 1.85 to 1.9 can be achieved. 

The present results suggest that impact volatilization can give rise to the 
greater density of the larger Jovian and Saturnian satellites as a result of vapori- 
zation of porous, water, ice, and possibly also NHg ice during accretion and all 
but the essentially silicate satellites can have accreted from the same primor- 
dial solar composition material. The final density of Jovian and Saturnian satel- 
lites then may just be the consequence of the amount of material in orbit 
around Jupiter and Saturn at a particular radius, and hence, ultimate satellite 
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4. SUMMARY AND CONCLUSIONS 

Shock wave data for polycrystalline, porous ice, and liquid water are avail- 
able over a sufficiently wide range of pressure. The peak shock state achieved 
upon accretionary impacts up to 2.6 km/sec, as well as, impacts of 20-30 
km/sec for cometary encounters are described using the impedance match 
method. A key calculation to determine the intersection of the pressure- 
volume-temperature-energy isentrope with the Kugoniot curve of water and ice 
so as to determine the entropy function along the Hugoniot. Knowledge of the 
entropy gain along the Hugoniot of ice determines whether upon rarefaction 
from a peak shock pressure will result in incipient melting, complete melting, or 
partial vaporization. The latter will occur upon passage of the isentropes 
through thermodynamic states near the vapor-liquid critical point (at 64 K, 221 
MPa and 3.1 m 3 /Mg). For ice the shock pressure required to induce incipient 
melting is 6 GPa, complete melting is 10 GPa, and passage through the vapor- 
liquid critical point is 22.5 GPa. These pressure values depend only slightly on 
the initial temperature of the water-ice. Upon isentropic release along the isen- 
trope passing through the critical point, the mean fraction of water vaporized is 
~0.5 which decreases to 0.4 upon cooling and expansion to below the vapor- 
liquid-ice I triple point. 

The vapor-liquid-ice I triple point at 61.08 kPa (6.1 mbar) and 273 K occurs 
at a point where the condensed phase density is 1 Mg/m 3 and the vapor phase 
density is 5 x 10 -a Mg/m 3 plays an important role in impact vaporization in -an 
environment with low partial pressure of H a 0. The minimum impact velocity 
onto an icy satellite surface required to induce the incipient melting is 2.1 and 
3.4 km/sec for silicate and icy projectiles, respectively. Complete melting of ice 
requires silicate or icy impactors traveling at ~3 and 4.4 km/sec, respectively. 
To obtain 0.4 to 0.5 mass fraction vapor upon passing through the liquid-vapor 
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critical point requires 4.5 to 7.2 km/sec impact velocities for silicate and ice 
projectiles, respectively. 

The Rice-Walsh-Bakanova equation of state in which internal energy or 
enthalpy is specified along a pressure-volume reference curve can be used to 
calculate other thermodynamic states such as along isentropes, at the same 
pressure, for different volumes and temperatures from the thermodynamic 
parameters (dH/ 3V)p or P(3V/3E)p. These parameters which are assumed to be 
only a function of P are fit to the singly or doubly shocked Hugoniot data for 
water, ice, and porous ice. 

The satellites of Jupiter and Saturn are classified into three groups accord- 
ing to density. The inner, silicate, satellites of Jupiter (Amalthea, Io, and 
Europe) the partially devolatilized (pd) satellites (Ganymede, Callisto, and 
Titan), and the smaller primordial (p) satellites (Mimas, Enseladus, Tethys, 
Dione, Rhea, and lapetus). The latter two groups have mean densities of 1.8 to 
1.9 Mg/m 3 and 1.0 to 1.4 Mg/m 3 . We infer that the densities of the P group also 
correspond to the densities of the particles in the rings of Saturn and Jupiter. 
We suggest that the p satellites are composed largely of primordial solar ratio 
material controlled by the oxygen/siiicon ratio (15.5) and the 
oxygen/magnesium ratio (17.4) of the primitive Jovian and Saturnian nebulae. 
The primordial solar ratios, if these correspond to objects having silicate and 
water ice compositions, would have densities of 1.1 to 1.2 Mg/m 3 . We have exam- 
ined the hypothesis that the large pd satellites achieved their higher densities 
(1.8 to 1.9 Mg/m 3 ) upon partial impact devolatilization during accretion starting 
with planetesimals of material similar in density, 1.27 Mg/m 3 , (60% ice, 40% sili- 
cate) to that of the smaller satellites. Accretion of planetesimals of such a com- 
position assuming impact volatilization at partial pressures of H 2 0 below 62 kPa 
in a porous regolith yields Ganymedean-sized planets with mean planetary 
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densities of only 1.5 Mg/m 3 . However, starting with planetesimais having densi- 
ties of 1.62 Mg/m 3 (43% H 2 0 and 57% silicate) give rise to mean satellite densi- 
ties of 1.9 Mg/m 3 for Ganymedean-sized objects. 

If impact volatilization of initially porous planetesimais is assumed our cal- 
culations demonstrate that starting with planetesimais of composition 54% H a 0 
and 46% silicate (mean density 1.35 Mg/m 3 ) and partial devolatilization during 
accretion occurs leads to objects having a Ganymedean radius of 2600 km and a 
mean density of 1.85 Mg/m 3 . 

The present results suggest a mechanism for the formation processes of the 
higher density partially devolatilized, satellites of Saturn and Jupiter. The accre- 
tion devolatilization process depends on the relative amount of material present 
in different zones orbiting in the nebular around each proto-planet and less on 
the assumed mean radial temperature at a given distance from the planet. 
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APPENDIX A. SHOCK TEMPERATURES AND THE WATER EQUATION OF STATE 

A stringent test of the complete thermal equation of state for fluids such as 
water are to compare predicted sound speeds and Hugoniot temperatures to 
measurements. Previously Bakanova et al. (1976) very favorably compared their 
predicted sound velocity along the water Hugoniot with measurements to ~60 
GPa. Shock temperatures were also calculated and compared with previous rela- 
tively low pressure data of Kormer (1968). Other work by Lyzenga and Ahrens 
(1980), Lyzenga et al. (1983a, b), Ahrens et al. (1982), Nellis et al. (1984), and 
Boslough et al. (1984) for a variety of solids and fluids demonstrate that upon 
comparing calculated and measured shock temperature Th versus shock pres- 
sure, P, the absolute value of Th at a given pressure was very sensitive to the 
total energy budget of the material (e.g. occurrence of phase changes) whereas 
a slope of the P versus Th relation is sensitive to specific heat. 

In Fig. 1A the available H 2 0 shock temperature data of Kormer (1968) and 
the more recent higher temperature and pressure data of Lyzenga et al. (1982) 
are shown relative to previous shock temperature calculations. As is evident 
from the figure, the calculation of Ree based on a repulsive parameter equation 
of state which does not take into account, the already mentioned, extensive 
pressure-induced ionization (Mitchell and Nellis, 1982) yields temperatures 
which are too high by ~600K, especially at higher pressures (above 50 GPa). 
Idee and Walsh’s calculations are too high largely because the data range they 
studied, up to ~25 GPa, did not encompass as much of the higher pressure range 
where subsequent measurements demonstrate that the Hugoniot curve is more 
compressive (especially above 30 GPa) (Mitchell and Nellis, 1982) than would be 
anticipated from simple extrapolation of the lower pressure data. Both the 
curve labeled Cowperthwaite and Shaw (1970) and a present calculations agree 
well with Kormer' s 1968 lower pressure results. Cowperthwaite and Shaw (1970) 
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use a thermal equation of state which was derived from a careful fit to the Lower 
pressure thermodynamic data, and took into account the strong temperature 
dependence of specific heat at constant volume. For the present calculation 
curve we employed the form of the Bakanova equation of state and used equa- 
tion 21 for specifying rj = P/ (dV/ 9E)p as a function of P. Instead of using a 
theoretical value of E i0 , T i0 , and V io to specify the foot of a series of isentropes at 
P i0 = 1 GPa in the low temperature and high pressure regime as do Bakanova et 
al. (1976), we used the experimentally based thermodynamic and equation of 
state data specified in Table 1A. Hence it is not surprising that our results agree 
so closely with those of Cowperthwaite and Shaw (1970). Finally we note that 
over the range from 65 to 80 GPa the Bakanova et al. (1976) equation of state is 
in excellent agreement with the recently measured shock temperatures. We 
thus conclude that this form of this equation of state provides an accurate 
description of the properties of shocked water and presumably ice up to 50 GPa 
when constrained by the thermodynamic data extending to ~1 GPa and ~1273K. 
At higher pressures > 50 GPa the theoretical description of the high tempera- 
ture - low pressure properties of water proffered by Bakanova et al. (1976) fits 
both the shock temperature data, the sound speed, and the reflected shock, and 
porous ice data. 
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Figure 1. Mean orbital radius versus satellite density for the Jovian and Saturni- 
an systems. 

Figure 2. Shock velocity versus particle velocity for ice at 253 and 263 K. 

Figure 3. Shock pressure versus particle velocity for water, ice, and porous ice. 
Curve for serpentinite, representative of planetary silicates after Marsh (1980). 
Escape velocity for satellites and planets indicated along the particle velocity 
axis. Shock pressures required for incipient melting (IM), complete melting 
(CM), and the critical point isentrope (CP) for ice are indicated. 

Figure 4. Bakanova parameter, versus, pressure. 

Figure 5. Principal Hugoniot and release isentropes for water in the pressure 
specific volume plane after Bakanova et al. (1976). 

Figure 6. Log 10 pressure (GPa) versus temperature phase diagram for H a 0. 
Hugoniot for ice and liquid water and release isentrope passing through liquid- 
vapor critical point (K) and water-vapor-ice I, triple point (CP) is indicated. 

Figure 7. Simplified Mollier enthalpy entropy diagram of water, after Bosnjako- 
vic, Renz, and Burow (1970). 

Figure 8. Hugoniot for ice centered at 70 and 263 K and release isentropes for 
incipient melting (IM), complete melting (CM), and critical point (CP) isentropes 
in pressure volume plane. Hugoniot of ice centered at 70 K inferred from ther- 
mal expansion data for ice collected in Fletcher (1970). 

Figure 9. Mass fraction of incident planetesimal material accreted, versus, 
planetary mass for pure H a 0 and 0.4 silicate + 0.6 H z 0 planet. A porous, ideal 
locking solid (similar to snow) is assumed on the surface of the growing planet. 
All the gravitational energy of accretion is deposited in this layer and gives rise 
to partial vaporization of ice at temperatures below 273.16 K and 61.03 kPa 
(6.108 mbar). 

Figure 10. Average uncompressed density of an icy planet versus planetary ra- 
dius. Light curves are for accretion of solid planetesimals and a regolith- 
covered planet. Heavy curve is for accretion of porous planetesimals and a solid 
planetary surface. 

Figure 1A. Shock temperature versus shock pressure. Theoretical calculations 
and experimental measurements. 
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Table 1. Characteristics of Major Jovian and Saturnian Satellites^ 


Satellite 

Mean Orbital Radius 
(10®km) 

Radius 

(km) 

Mass 

(lO^kg) 

Mean Density 
(Mg/m 3 ) 

Satellite 

Class. 

J-5 Amalthea 

1.8 

(96) 

(0.13) 

(3.5) 

(silicate) 

Jl-Io 

4.2 

1815±10 

892 

3.55 

silicate 

J2-Europa 

6.7 

1569±10 

487 

3.04 

silicate 

J3-Ganymede 

11 

2631±10 

1490 

1.93 

p.d. 

J4-Callisto 

19 

2400±10 

1075 

1.83 

p.d. 

Si-Mi mas 

1.9 

195±5 

0.37 

1.44±0. 18 

P 

S2-Enceladus 

2.4 

250±10 

0.72 

1. 16±0.55 

P 

S3-Tethys 

2.9 

525±10 

6.06 

1.2l±0. 16 

P 

S4-Dione 

3.8 

560±10 

10.30 

1.43±0.06 

P 

S5-Rhea 

5.3 

765±10 

24.38 

1.33±0.09 

P 

S 8- Titan 

1Z2 

2560±26 

1335.25 

1.88±0.01 

p.d. 

S8-Iapetus 

35.6 

720±20 

18.78 

1.16±.09 

P 


M Data compiled from Morrison (1982), Stone and Miner (1982), Lunine and Stevenson (1982). 



Table 2. Shock-Particle Velocity Data for Water, Ice, and Porous Ice 


Initial Material- 

Initial Density 

C 9 (km/sec) 


Aa 

Range of 

Ref. 


(Mg/m 3 ) 



(sec /km) 

u (km/ sec) 


Water (294-298K) 

0.9979±0.0004 

2.393 

1.333 

- 

1.5<u<7.1 

(a) 

Water (293) 

0.99823^ 

1.50 

2.0 

-0.107 

0<u<4.0 

(c) 

Water (293) 

0.99 823 W 

3.2 

1.144 

- 

u>4.0 

(c) 

Ice (258) 

0,915 

1.57 

1.485 

- 


(c) 

Ice (283) 

0.92S0W 

1.317 

1.526 

- 


(d) 

Porous Ice (258) 

0.80 

0.74 

1.425 

- 


(o) 

Porous Ice (258) 

0.35 

0.0 

1.425 

- 


(c) 


W Mitchell and Nellis (1982) 

ft) Nominal value Weast (1982) 

ft) Bakanova et al. (1978) i 

ft) present fit to 18 data points (8 data, Bakanova et al., 1976; 8 data, Anderson, 1988; 1 point, Larsen 
et al., 1972; 1 point, Ahrens and Gaffney, 1978). 



Table 3. Relative cratering rates on the Galilean satellites for mean 
encounter velocity of 8 km/sec (Smith et al., 19?9b) 


Satellite 

Flux concentration 
factor by Jupiter gravity 

Mean Impact 
Velocity (km/sec) 

Relative Crater 
Production Rate 

Callisto 

3.3 

14 

1 

Ganymede 

5.0 

18 

2.0 

Europa 

7.3 

21 

2.8 

Io 

11.1 

26 

5.3 



Table 4. Critical Isentropes in the Water System 


Critical Isentrope 

Determining 

Entropy* 


State 

(kJ/kgK) 

Incipient Melting 

Triple Point 

2.21 

(IM) 

Ice 111-Ice V -Liquid 


Complete Melting 

Triple Point 

3.16 

(CM) 

Ice Hoe Ill-Liquid 


Liquid-Vapor 
Critical Point (CP) 

Liquid-Vapor 

7.95 

1 GPa and 1273K 

none 

6.76 


State at 1 GPa 

Phase 

Temperature 

(K) 

Enthalpy* 

(kJ/kg) 

Volume 

(Mg/m a ) 

Internal Energy 
(kJ/kg) 

Ice VI 

400 

1140 

0.735 

405 

Ice VI+ 
Liquid 

400 

1400 

0.736 

614 

Liquid 

1017 

40BD 

1.085 

2995 

Vapor 

1273 

5022 

1.2532 

3769 


* referenced to 1 bar and OK 


Table 5. Critical Release States, Ice Hugoniots 



70K 

263K 

Isentrope Entropy W 

(kJVkgK) 

T(K) 

P(GPa) 

V(m 3 /Mg) 

TOO 

P(GPa) 

V(m 3 /Mg) 

01 

2.21 

378 

7.6 

0.594 

363 

6.2 

0.616 

CM 

3.16 

410 

10.8 

0.562 

298 

9.6 

0.576 

CP 

7.95 

1756 

25.5 

0.509 

1779 

26.7 

0.500 


(a) referenced to 1 bar and OK 



Table 1A. Shock Temperature, Liquid Water 


Isentropic 

Centered State, 

1 GPa 

Kugoniot State 

E(kJ/kg 

V(m 3 /Mg) 

T(K) 

E(kJ/kg) 

V(rr.VMg) 

T(K) 

P(GPa) 

3166W 

1.1 1ft) 

1073ft) 

8766 

0.484 

1964 

31.5 

3769ft) 

1.25ft) 

1273< a ) 

10,987 

0.460 

2447 

37.9 

(c) 

0.868ft) 

428ft) 

(c) 

0.595 

584 

10.9 

(c) 

1.014ft) 

1236ft) 

(c) 

0.443 

2555 

44 

(c) 

1.19ft) 

2143ft) 

(c) 

0.382 

4703 

75 

(c) 

1.381ft) 

3000ft) 

(c) 

0.344 

6959 

107.4 


(a) Bosnjakovic et al., 1970 

(b) from Burnham et al,, 1969 

(c) not calculated 

(d) Bakanova et al., 1976 


